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Abstract: In many situations, in all branehes of physios, one enoounters power-like 

behavior of some variables whieh are best desoribed by a Tsallis distribution eharaoterized 
by a nonextensivity parameter q and soale parameter T. However, there exist experimental 
results whieh can be described only by a Tsallis distributions which are additionally 
decorated by some log-periodic oscillating factor. We argue that such a factor can originate 
from allowing for a oomplex nonextensivity parameter q. The possible information conveyed 
by sueh an approaeh (like the oeeurrenee of eomplex heat eapaeity, the notion of eomplex 
probability or oomplex multiplieative noise) will also be disoussed. 

Keywords: soale invariance, log-periodio oscillation, complex nonextensivity parameter, 
complex multiplicative noise 


1. Introduction 


In many situations, in all branches of physics, one encounters behavior of some variables X which 
become pure power distributions for large values of X and exponential for X —0. Because of this they 
are known as power-like distributions and in many cases they are identified with a Tsallis distribution 
[ 1 ], 


1 - (1 - g) 




n i/(i-9) 


F{X) = A 


( 1 ) 
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characterized by a seale parameter T and parameter q known as nonextensivity parameter (^4 is 
normalization) k Obviously, for X —)■ 0 distribution (1) beeomes the usual Boltzmann-Gibbs exponential 
formula with temperature T, but it beeomes pure exponential (i.e., BG) also for g ^ 1. For q ^ 1 and 
large values of X it beeomes pure power distribution not sensitive to seale parameter T. 

To fully reeognize the nontrivial eharaeter of distribution (1), one must realize that, usually, in 
different parts of phase spaee of the variable X, one eneounters (or, rather, one expeets) a dominanee 
of different (if not eompletely disparate) dynamieal faetors. This is best seen in the proeesses of 
multipartiele produetion at high energies (the best known to us). They will serve here to exemplify 
our further eonsideration eoneeming some speeifie log-periodie oseillations, apparently visible in sueh 
proeesses, whieh must be therefore somehow hidden in the original distribution (1). 

Before proeeeding further, we shall briefly summarize the present status of applieation of Tsallis 
distributions in this eontext, eoneentrating only on multipartiele produetion proeesses. They eomprise 
of many different meehanisms in different parts of phase spaee. Limiting ourselves only to partiele 
produetion in the eentral rapidity region and to distribution of their transverse momenta pt, it is 
eustomary to divide this produetion into independent soft and hard proeesses populating different parts 
of the transverse momentum spaee^ separated by a momentum seale po. As a rule of thumb, the speetra 
of the soft proeesses in the low-pr region are (almost) exponential, F(pT)~exp(—p-^/T), and are 
usually assoeiated with the thermodynamieal deseription of the hadronizing system. The px speetra 
of the hard proeess in the high-p^ region are regarded as essentially power-like, F{pt)^Pt^, and are 
usually assoeiated with the hard seattering proeess (for relevant literature eoneeming both parts see [2]). 
However, it was very soon reeognized that both deseriptions eould be replaeed by a simple interpolating 
formula [3], 

F{pt) = a(i + ^) , ( 2 ) 

V PoJ 

that beeomes power-like for high px and exponential-like for low px- The reasoning was that for high px, 
where we are usually negleeting the eonstant term, the seale parameter po beeomes irrelevant, whereas 
for low Px it beeomes, together with the power index n, an effeetive temperature T = po/n. The same 
formula re-emerged later to beeome known as the QCD-based Hagedorn formula [4]. It was used for the 
first time in [5] and beeame one of the standard phenomenologieal formulas for p^ data analysis [6-10]. 
In the mean time it was realized that both formulas are, in faet, identieal onee 

n = - and po = nT, (3) 

q - 1 


* The reason being the fact that Eq. (1) is also emerging from nonextensive statistical mechanics [1]. 

^ A few words of definition concerning this phase space is necessary. A produced particle has some momentum p = 
[pl,Pt]- Its longitudinal part, p^, is defined as parallel to the axis of collision, its transverse part, px as perpendicular to 
that axis. They are defined by means of rapidity y variable, y = ^ In -§ 3 ^, as, respectively, p = \p\ = m sinh y, whereas 
energy of particle, E = ^/rnXf-pA = m coshy. Central rapidity means y = 0. In what follows, our X from Eq. (1) will 
be identified with transverse momentum, X = px- 
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and therefore they ean be used interehangeably^ 

This distribution is usually used in a thermodynamieal eontent in whieh the seale parameter T is 
identified with the usual temperature (although such identification cannot be solid [25]) and with a real 
power index n = l/(g — 1) (or real nonextensivity parameter q). Actually, a Tsallis distribution can 
be regarded as a generalization to real power n (or q) of such well known distributions as the Snedecor 
distribution (with n = (z/ + 2)/2 and integer u, which for —)■ cxo it becomes exponential distribution). 

In [26] we investigate the case when g is a complex number. We shall review our results in this field 
in the next section adding examples where log-periodic oscillations occur at different energies and for 
different collision systems. In Section 3 we discuss the possible consequences of complex nonextensivity 
parameter including some new recent developments in this field (as complex probability and complex 
multiplicative noise). The final section contains our conclusions and summary. 

2. Log-periodic oscillations in Tsallis distribution - complex power index 

Recently, the experiments [8-10] at the Large Hadron Collider (LHC) at CERN provided new data in 
a very large domain of transverse momenta, pt, phase space. They turned out to be extremely interesting 
because of the following: 

• They allow us to test the standard Tsallis formula, Eq. (1), over ~ 14 orders of magnitude. As 
can be seen in Eig. la, the observed pt distributions of secondaries produced in proton-proton 
collisions in these experiments can be very well reproduced (cf. also [23]^. 

• And, what is of special importance to us, they disclose some features which suggest a departures 
from the single form of Eq. (1), cf.Eigs. 1 b-c.. Apparently they could not be seen in previous 
experiments because they seem to be connected with rather large values of transverse momenta, 
not available earlier.. 

However, whereas fits to Eq. (1) look pretty good, closer inspection shows that the ratio of data/fit is not 
flat. It shows some kind of visible oscillations, cf. Eig. lb. These are the oscillations we have mentioned 
before. 

It turns out that these oscillations cannot be compensated, or erased, by any reasonable change of 
fitting parameters. Moreover, they are visible by all three experiments CMS, ATEAS, AEICE. The only 
condition for such an effect to be visible is that the experiment covers a sufficiently large domain of 
transverse momenta pt, cf. Eig. lb. It is also seen at all energies covered by these experiments, cf. Eig. 
Ic. And, finally, as Eig. Id shows, this effect is also visible (and is even more pronounced) in nuclear 
collisions. When taken seriously, it turns out that to account for these oscillations one has to "decorate" 


^ Both Eqs. (1)) and (2) have been widely used in the phenomenological analysis of multiparticle productions, including 
situations where the nowadays observed spectra extend over many orders of magnitude, [6-23]. Up to now such possibility 
of testing Tsallis distribution offered only cosmic ray fluxes, cf. [24]. 

^ These secondaries were produced at midrapidity, i.e., for y ~ 0 for which, for large transverse momentum, pT > m 
(where m is the mass of the particle), one has that, approximately, the energy of particle, E ~ pr, i.e., it practically 
coincides with pT- 
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distribution / {pt) from Eq. (1) (i.e., one has to multiply it) with some log-periodie oseillating faetor. It 
is is usually taken in the form [29]: 

R{E) = a + 6cos [c\n{E + d) + f]. (4) 


Figure 1. (Color online) Examples of log-periodic oscillations, (a) dN/dpT for highest 
energy 7 TeV, the Tsallis behavior is evident. Only CMS data are shown [8], others behave 
essentially in an identical manner. (6) Eog-periodic oscillations showing up in different 
experimental data like CMS [8] or ATEAS [9] taken at 7 TeV. (c) Results from CMS [8] for 
different energies, (d) Results for different systems (p + p collisions compared with Pb + Pb 
taken for 5 % centrality [27]. Results from AEICE [28] are very similar. Eits for p + p 
collision at 7, 2.76 and 0.9 TeV are performed with g = 1.139 + i-0.0385, 1.134 + i-0.0269 
and 1.117 + i ■ 0.0307, respectively. Eit for central Pb + Pb collisions at 2.76 TeV is done 
with q = 1.135 + i ■ 0.0321. See text for more details. 






Before proceeding any further let us remember that such log-periodic oscillations are widely know 
in all situations in which one encounters power distributions. In fact, such behavior has been found 
in earthquakes [30], escape probabilities in chaotic maps close to crisis [31], biased diffusion of 
tracers on random systems [32], kinetic and dynamic processes on random quenched and fractal 
media [33], when considering the specific heat associated with self-similar [34] or fractal spectra [35], 
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diffusion-limited-aggregate elusters [36], growth models [37], or stoek markets near finaneial erashes 
[38], to name only a few examples. However, in all these eases the basie distributions were seale free 
power law, without any seale parameter (here T) and without a eonstant term governing their X < nT 
behavior. 

In the eontext of nonextensive statistieal meehanies log-periodie oseillations have first been observed 
and diseussed while analyzing the eonvergenee dynamies of 2 ;-logistie maps [39]. In this paper we 
shall propose another way of introdueing sueh oseillations to Tsallis distributions. It will be based on 
allowing the power index n (or nonextensivity parameter q) in a Tsallis distribution to beeome eomplex. 
For eompleteness of the presentation we start from the simple pure power law distribution, 

0{x)=C- x-^. (5) 


This funetion is seale invariant, i.e., 

O (Ax) = pO (x), (6) 

with m = — In yu/ In A. However, beeause 1 = exp {i^nk), one ean as well write that 

/xA™' = 1 = exp (z27r/c), /c = 0,1,.... (7) 


It means therefore that, in general, the index m ean beeome eomplex. 


m = 


In /i 

MA 


-f ^ 


2'Kk 

hTI' 


( 8 ) 


As will be obvious from further, general eonsiderations, sueh a form of the power index results in R as 
given by Eq. (4) when one only keeps /c = 0,1 terms (whieh is the usual assumption eustomary applied 
in all applieations [29-33]). 

However, Tsallis distribution is only a power-like, not a power distribution. Therefore, to explain the 
origin of sueh a dressing faetor in this ease one has to find a right variable in whieh the real sealing holds. 
We start from the observation that, whereas the Boltzmann-Gibbs (BG) distribution. 



(9) 


eomes from the simple equation. 


df{E) 

dE 



( 10 ) 


with the seale parameter T being eonstant, the same equation, but with variable seale parameter in the 
form 


T = T{E) = To + -, 
n 


( 11 ) 


(known as preferential attachment in networks [14,40]^ ), 


df{E) 

dE 


1 


f\E) 


1 

Tq -f E / n 


fiE), 


( 12 ) 


5 


It is worth recalling here that this very same form, T{E) = Tq -f (1 — q)E, also appears in [22] within a Fokker-Planck 
dynamics applied to the thermalization of quarks in a quark-gluon plasma by collision processes. 









Entropy 2014 , xx 


6 


results in the Tsallis distribution 




riTo \ nTo 

We shall write now Eq. (12) in finite differenee form, 


fiE + 5E) = 


( 13 ) 


nT + E 


(14) 


In praetieal sense this means a first-order Taylor expansion for small 5E « E (from Eq. (14) on, 
we use T instead of Tq). We shall now consider a situation in which 5E always remains finite (albeit, 
depending on the value of the new scale parameter a, it can be very small) and equal to 


5E = anT{E) = a{nT + E). 


(15) 


Because one expects that changes 5E are of the order of the temperature T, the scale parameter must be 
limited by 1/n, i.e., a < 1/n. In this case, substituting (15) into (14), we have. 


f[E + a{nT + E)] = (1 - an)f{E). 

Expressing Eq. (16) in a new variable x. 


X =1 + 


E 


(16) 


(17) 


we recognize that the argument of the function on the left-hand side of equality (16) is 

E + a{nT -I- £■) = (1 -f a)xnT — riT, 
while the argument of the function on its right-hand side is 


E = xnT — nT. 


Notice that, in comparison with the right-hand side, the variable x on the left-hand side is multiplied by 
the additional factor (1 -I- a). This means that, formally, Eq.(16), when expressed in x, corresponds to 
the following scale invariant relation: 


g[{l + a)x] = {l-an)g{x). (18) 

This means than that, following the discussion after Eq. (6), its general solution is a power law, 

g{x)=x-^\ (19) 

with exponent irik depending on a and acquiring an imaginary part, 

ln(l — an) 271 

' ik- 


ln(l -f a) ln(l -f a) 


mk = 


( 20 ) 
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The special case of A; = 0, i.e., the usual real power law solution with mo corresponding to fully 
continuous scale invariance®, recovers in the limit a —)■ 0 the power n in the usual Tsallis distribution. 
In general one has 


g{x) = ^ Wk ■ Re {x = x ^ wu ■ cos [Im {mk) ln(x)]. (21) 

fc =0 fc =0 

One therefore obtains a Tsallis distribution decorated by a weighted sum of log-oscillating factors (where 
X is given by Eq. (17)). Because usually in practice we do not a priori know the details of the dynamics 
of processes under consideration (i.e., we do not known the weights Wk), for fitting purposes one usually 
uses only fc = 0 and A; = 1. In this case one has, approximately. 


g{E) ~ 


E 


Wq + Wi cos 


27r 


ln(l + a) 


In 


E \ 


( 22 ) 


and reproduces the general form of a dressing factor given by Eq. (4) and often used in the literature 
[29]. In this approximation the parameters a, b, c, d and / from Eq. (4) get the following meaning: 

a wq 271 , „ „ 271 

b wi ’ ln(l + a) ’ 


d = nT, f = 


ln(l + a) 


In(nT). 


(23) 


In fact this is not the most general result for in our derivation, Eqs.(15)-(18)), we have so far only 
accounted for a single step evolution. In real situation one should expect to have a whole hierarchy of 
evolutions. In such a case consecutive steps of evolution are connected by: 


Ei — Ei_i + CTj-i {nT + Ei_i ), 


(24) 


each with its own scale parameter ctj. In the simplest situation, neglecting any fluctuations of consecutive 
scaling parameters, i.e., assuming that all Oj = a, one has that after n steps 

nT + E/fj = (1 + o;)^ [nT + Eq) . (25) 


This means that, in general, Eq. (18) should be replaced by a new scale invariant equation: 

g [(1 + a)'"x] = (1 - an)'^g{x). 


(26) 


Whereas this equation does not change the slope parameter tuq, it significantly influences the frequency 
of oscillations which are now k times smaller. 


27r 

fi:ln(l + a) 


(27) 


(in Eq.(26) A = (1 + a)^ and p = [1 — an)'^; the slope parameter mo = — In /r/ In A is independent of k, 
whereas the frequency of oscillations, 27t/ In A, decreases with k as l/n). Eor more complex behavior of 
intermediate scale parameters a* one gets more complicated expressions (we shall not discuss this here). 


In this case power law exponent toq still depends on a and increases with it roughly as mo — n + j(n + l)a -\- 
^ (4n^ + 3n — l) + -^ (6n^ + 4n^ — n -I- l) -I- -Notice also that a < 1/n. 


6 
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3. Other consequences of complex nonextensivity parameter 

There are other eonsequenees of allowing the parameter m to be complex. In what follows we shall 
discuss shortly three examples: complex heat capacity , complex probability and complex multiplicative 
noise. 


3.1. Complex heat capacity 


The complex power exponent in the Tsallis distribution, m = m' + i ■ m", means that 


q-l = 


1 

m 


m ' 
mP 


+ i 



(28) 


As shown in [18] (cf. also [14,15,41]), the nonextensivity parameter q can be treated as a measure of 
the thermal bath heat capacity C with 

C = - = m! + im". (29) 

g- 1 

The complex nonextensive parameter q must therefore have some profound consequences because now 
the corresponding heat capacity becomes complex as well. As a matter of fact, such complex (frequency 
dependent) heat capacities (or generalized calorimetric susceptibilities) are known in the literature [45] 
and are usually written in the form 

C = Coo+ - iu:t). (30) 

1 + [ojtY 


Here Coo is the heat capacity related to the infinitely fast degrees of freedom of the system as compared 
to the frequency cu, and Cq is the total contribution at equilibrium (the frequency is set to zero) of the 
degrees of freedom, fast and slow, of the sample. The time constant r is the kinetic relaxation time 
constant of a certain internal degree of freedom. 

These complex heat capacities are known as dynamic heat capacities and are intensively explored 
from both experimental and theoretical perspectives. It is expected that dynamic calorimetry can provide 
an insight into the energy landscape dynamics, cf., for example, [46-49]. Usually one associates the 
imaginary part of linear susceptibility with the absorption of energy by the sample from the applied 
field. 

In the case of temperature fluctuations 5T{t) the deviation of the energy from its equilibrium value 
5U{t) is, for a certain linear operator C{t), some linear function of the corresponding variation of the 
temperature, 

5U{t) = C5T{t). (31) 


If the temperature of the reservoir changes infinitely slowly in time, then the system can keep up with 
any changes in the reservoir and its susceptibility is just the specific heat of the system Cy- However, 
in general, the behavior of the system is described by a generalized susceptibility Cv{oj), which can be 
called the complex and u-dependent heat capacity of the system. The change in the energy of a system 
in the field of the thermal force can be represented by 


5U{t) = 


L {t') 5T {t - t') dt', 


(32) 






Entropy 2014 , xx 


9 


where L (t') is the response function of the system describing its relaxation properties given by <f)(t) = 
L {t') dt'. Taking the Fourier transform one gets 


SU{uj) = Cv{uj)ST{uj), 


(33) 


where 

Cvico) = j (34) 

is the generalized susceptibility of the system and is called the complex heat capacity. In practice, the 
frequency dependent heat capacity is a linear susceptibility describing the response of the system to the 
small thermal perturbation (occurring on the time scale l/ce) that takes the system slightly away from 
the equilibrium . 

A complex Cy(a;) means that 5U and 5T are shifted in phase and that the entropy production in 
the system differs from zero [49]. The corresponding fluctuation-dissipation theorem for the frequency 
dependent heat capacity was established in [48]. According to this result, the frequency-dependent heat 
capacity may be expressed within the linear response approximation as a linear susceptibility describing 
the response of the system to arbitrarily small temperature perturbations away from equilibrium, 

Cv(u) = j" dte-<“‘{U(,0)U(t)) (35) 


(the u denotes frequency with which temperature field is varying with time). 

The above results for heat capacity can now be used to a new phenomenological interpretation of the 
complex q parameter discussed before. Namely, one can argue that 


Var{T) . S{T) 

(T)2 


(36) 


were 

S{T) = uj J Cov[T{0),T{t)]e-^‘^^dt (37) 

is the spectral density of temperature fluctuations (i.e., the Fourier transform of the covariance function 
averaging over the nonequilibrium density matrix). 

We would like to stress at this point that, in a sense, Eq. (36) can be regarded as a generalization of our 
old proposition for interpreting g as a measure of nonstatistical intrinsic fluctuations in the system [43,44] 
(which corresponds to the real part of (36)) by adding the effect of spectral density of such fluctuations 
(via the imaginary part of (36)). Notice that (36) follows from (29) and the relation U = CyT, allowing 
to write (35) in the form of (36). 


3.2. Complex probability 

From the point of view of superstatistics [42,43], in our particular case complex parameter q 
corresponds to a complex probability distribution. Namely, one uses the property that gamma-like 
fluctuation of the scale parameter T in an exponential BG distribution (9) results in the g-exponential 
Tsallis distribution (1) with g > 1. The parameter g is given here by the strength of these fluctuations, 
g = l + Var{X)/ < X > 2 . From the thermal perspective, it corresponds to situation in which the heath 
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bath is not homogeneous but has different temperatures in different parts, which are fluctuating around 
some mean temperature Tq. It must be therefore described by two parameters: a mean temperature Tq 
and the mean strength of fluctuations given by q. 

We now perform the same procedure, but using two gamma distributions, one with a real power index, 
mo — 1, and one with a complex power index, mo + imi — 1, 

(" 1 ) “P (-"!) + 

mo+imi —1 / 'X' \ 

exp ■ ( 38 ) 

As the result one gets a complex distribution (complex pdf): 

/ oo / F \ / F \ 

f(E)g(l/T)d(l/T) = Cwo + j ’ 

the real part of which is pdf in form of a Tsallis distribution decorated with log-periodic oscillations of 
the type of Eq. (22), 


+tei 


r (mo -f imi) 


T I 

niQ n— 


Re[hq{E)\ = C (^ +-iwo + wicos -f 

The complex pdf has a number of interesting properties [50,51]. It plays an important role in the 
interference among resonance states during scattering experiments. It is associated with the phase 
of the resonance channel probability amplitudes (in non-Hermitian quantum mechanics). In wireless 
communication systems it is generated by a superposition of finite random variables and usually involves 
the movement, scattering, diffusion or diffraction. The imaginary part is proportional to the degree of 
the correlation. The imaginary part is then a function of a correlation coefficient or other parameters 
that state the degree of the relationship of each individual random variable of the superposition of the 
random variable having a complex pdf. The real and imaginary part have diverse properties, i.e. one for 
real valued pdf and the other for elementary correlation, respectively. 

It is interesting to note that entropy 



H 


'a In a + i ■ 6 In h)dxidx 2 


(41) 


corresponding to complex joint probability. 


/ (xi, X2) = a (xi, X2) -t- i ■ 6 (xi, X2), 


(42) 


consists of two components: 


H, = 


alnadxi(ix2, H2 = — b\nbdxidx2', H = \Hi + iH2\y Hf + H 2 > Hi. (43) 


The imaginary part of entropy is proportional to the degree of incompatibility of the correlated stochastic 
processes. The incompatibility increases the entropy of correlated stochastic processes. 
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3.3. Complex multiplicative noise 


It is known that multiplicative noise leads to a Tsallis distribution [44]. It is then natural to expeet 
that multiplieative eomplex noise should result in eomplex q and in log-periodie oseillations in Tsallis 
distributions. It ean be defined by a Langevin equation 

(I'D 

-E + = ^(t), where 7 (f) = 70 (f) + iyi. (44) 

The resulting distribution [44] is now 


fip) 



where T 


2Var{0 

(7) 


Q = 


2Var(y) 

(7) ■ 


(45) 


The parameter q is now eomplex beeause (7) is eomplex. Even more importantly, (g — 1)/T = 
Var{y)/Var{^) is real (it tends to zero for q —)■ 1). This is beeause the eomplex term 71 added to the 
noise is eonstant. Notiee that we eould just as well replaee in Eq. (45) (q — 1) {p^/T) by (p^/po) where 
Po = Var{^)/Var{y). The examples and diseussion of the systems eharaeterized by the appearanee of 
"imaginary" multiplieative noise terms in an effeetive Eangevin-type deseription ean be found in [52]^. 


4. Summary and conclusions 

In may places in physics, and especially in the realm of high energy multiparticle production processes 
we are particularly interested in, it became a standard procedure to fit the data on transverse momentum 
distributions by means of the quasi-power Tsallis formula. The usual interpretation in such cases is 
that the scale parameter T is a kind of "temperature" whereas additional nonextensivity parameter q is 
describes intrinsic, nonstatistical fluctuations existing in the system [1 1-22,24,42-44,53]. However, with 
increasing range of transverse momenta measured in recent experiments [ 8 - 10 ] two things happened: 

(f) That they still can be fitted by the same formula (which came as surprise because fits now cover 
~ 14 orders of magnitude of the measured cross sections [2,23]). 

{a) That new data revealed weak but persistent oscillation of log-periodic character (discussed already 
shortly in [26]). 

If taken seriously, such log-periodic structures in the data indicate that the system and/or the underlying 
physical mechanisms have characteristic scale invariant behavior. This is interesting as it provides 
important constraints on the underlying physics. The presence of log-periodic features signals the 
existence of important physical structures hidden in the fully scale invariant description. It is important 
to recognize that Eq. (12) represents an averaging over highly ’non-smooth’ processes and, in its present 
form, suggests rather smooth behavior. In reality, there is a discrete time evolution for the number of 
steps. To account for this fact, one replaces a differential Eq. (10) by a difference quotient and expresses 


^ In fact, this is not exactly Tsallis formula from Eq. (1). To get it one has to allow for correlation between noises and drift 
term due to additive noise, i.e., for Cov{C l) 0 and (^) 7 ^ 0 (see [53] for details). One obtains then Eq. (1) but with, in 
general, complex T = T{q). We shall not discuss it here. 
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dt as a discrete step approximation given by Eq. (15) with parameter a being a charaeteristie seale 
ratio. It can also be shown that diserete seale invarianee and its associated eomplex exponents ean appear 
spontaneously, without a pre-existing hierarehical strueture. Finally, a complex nonextensivity parameter 
promises new perspeetives in future phenomenologieal applieations being eonneeted to eomplex heat 
eapaeity, to notion of eomplex probability or to eomplex multiplieative noise, to mention only a few 
examples diseussed shortly in our paper. 
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